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SUMMARY 


A.  TECHNICAL  PROBLEM 

This  project  is  one  part  of  a  larger  effort  directed  at  extending  packet  switching 
technology  into  the  area  of  radio  communications.  The  technical  investigations  are 
focused  on  a  Packet  Radio  Network  to  serve  fixed  and  mobile  digital  terminals  in  a 
tactical  command  and  control  environment.  The  objectives  of  the  Collins  investigations 
are: 


•  Perform  research  covering  the  application  of  radio  frequency  technology  to 
packet  switched  communications. 

•  Participate  in  the  overall  ARPA  Packet  Communications  Technology  Program 
under  the  guidance  of  the  Packet  Radio  Communications  Working  Group. 

•  Develop  experimental  equipments  to  support  the  Packet  Radio  System/Network 
experiments 

•  Recommend  a  system  architecture  to  provide  terminal-to-terminal  security 

|  for  the  packet  switched  radio  network. 

B.  GENERAL  METHODOLOGY 

The  approach  to  the  packet  radio  investigations  is  a  combination  of  theoretical 
studies,  equipment  design  efforts,  and  laboratory  experiments.  Work  has  been 
directed  at  each  of  the  following  areas: 

•  Experimental  Equipment  Design  and  Build.  Detailed  design  and  assembly  of  the 
experimental  repeater  described  in  the  last  quarterly  report,  ''Technical  Plan 

1  for  the  Packet  Radio  Experimental  Repeater,  "  July  1974. 

•  Modulation  and  Detection.  Investigation  of  modulation  schemes  and  signal 
processing  techniques  as  related  to  contention,  channelization,  and  coexistence 
properties. 

•  Communications  Security.  Evaluation  of  the  impact  of  encryption  and  traffic 
flow  security  requirements  and  disciplines  on  system  design. 

The  above  studies  and  design  efforts  are  organized  as  an  interactive  process  to 
achieve  a  proper  balance  between  system  performance  and  the  constraints  of 
available  technology. 

I 


iv 


C.  TECHNICAL  RESULTS 

This  quarterly  report  treats  three  primary  topics  which  have  been  investigated 
during  the  reporting  period  in  support  of  the  experimental  packet  radio  repeater 
development.  Specific  technical  results  reported  include  the  following: 

•  Pseudonoise  Spread  Spectrum  Modulation.  Detailed  consideration  of  pseudonoise 
spread  spectrum  (PNSS)  techniques,  including  coexistence  considerations. 

•  Code  Selection.  Description  of  an  optimization  program  for  selecting  the 
Golay  code  sequences  to  be  used  for  MSK  modulated  code  spread  spectrum 
signaling  in  the  experimental  repeater. 

•  Capture  Model  Description.  Characterization  of  the  ability  of  the  experimental 
repeater  receiver  design  to  distinguish  and  lock  onto  one  signal  from  an  ensemble 
of  several  simultaneous  signals. 

In  addition,  substantial  progress  has  been  made  in  the  detailed  design  and  development 
of  the  experimental  packet  radio  repeater.  Breadboards  of  most  circuit  functions 
have  been  built  and  are  being  tested.  Printed  circuit  boards/modules  are  under 
development. 

The  report  of  work  done  in  communications  security  is  reported  in  other  documents. 

D.  PLAN  FOR  NEXT  QUARTER  ACTIVITIES 

The  emphasis  during  the  next  quarter  will  be  on  the  continued  design  and  integration 
of  the  experimental  repeater.  Circuit  design  and  layout  will  be  completed;  modules 
and  circuit  boards  will  continue  to  be  fabricated,  assembled,  and  tested;  and  the 
radio  and  digital  subsections  of  the  experimental  repeater  will  be  integrated.  Extensive 
tests  will  be  run  at  this  level  to  thoroughly  evaluate  the  overall  performance  of  the 
experimental  repeater  design  prior  to  final  assembly,  test,  and  integration  of  the 
initial  experimental  units  into  the  packet  radio  experimental  test  bed. 
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Section  1 


Introduction 


A  spread  spectrum  mode  of  signaling  has  been  selected  for  the  experimental  repeater  for 
the  ARPA  packet  radio  project.  The  choice  was  made  because  the  processing  gain  associated 
with  the  expanded  bandwidth  resulting  from  the  use  of  spread  spectrum  techniques  can  be 
used  to  provide  a  degree  of  "capture"  for  the  multiple  access  channel  or  some  protection 
against  interfering  signals.  The  use  of  the  spread  spectrum  techniques  also  reduces  the 
effect  of  the  packet  radio  signals  on  other  systems  operating  in  the  same  general  frequency 

band. 

This  report  presents  the  results  of  several  studies  of  pseudonoise  spread  spectrum  modu¬ 
lation  techniques  and  their  use  in  the  experimental  repeater  now  under  development  for  the 
ARPA  packet  radio  project  Section  2  is  a  detailed  discussion  of  pseudonoise  spread 
spectrum  modulation.  Performance  advantages  with  respect  to  operation  in  the  presence 
of  strong  interfering  signals,  such  as  jamming  signals,  coexistence,  multipath  rejection, 
and  multiple  access  to  a  single  channel,  are  discussed.  Section  3  is  a  description  of  an 
optimization  program  developed  for  selecting  Golay  code  sequences  to  be  used  for  MSK 
modulated  code  spread  spectrum  signaling  in  the  experimental  repeater.  Finally,  Section  4 
contains  a  characterization  of  the  expected  signal  capture  capability  of  the  experimental 
repeater,  that  is,  the  ability  of  the  receiver  to  distinguish  and  lock  onto  one  spread  spectrum 
radio  signal  from  an  ensemble  of  several  simultaneous  signals. 

In  addition  to  the  work  described  above,  s:  bstantial  progress  has  been  made  on  the  develop¬ 
ment  of  the  experimented  packet  radio  repeater.  Breadboards  of  most  of  the  radio  section 
circuits  have  been  completed  pursuant  to  the  technical  plan  (published  in  the  last  quarterly 
report)  and  are  undergoing  test.  Completed  breadboards  include  the  encoder/modulator,  the 
rf  head  (up-converter,  power  amplifier,  T/R  switch,  and  down-converter),  the  noncoherent 
age  the  coherent  age,  EOP/data  detector,  and  the  bit  sync.  In  addition,  printed  circuit 
board  modules  are  under  development  for  most  of  these  functions.  A  completed  prototype 
module  is  available  for  the  frequency  generator.  With  respect  to  the  digital  section,  printed 
circuit  boards  for  the  random  access  memory  and  the  programmable  read  only  memory 
have  been  built  and  tested. 

The  following  modules  have  been  designed  and  are  in  printed  circuit  board  layout:  CPU 
(including  timing),  address  decode,  DMA  channels,  and  the  receive  interface.  The  repeater 
software  development  is  proceeding  according  to  plan.  The  following  program  elements 
are  designed  and  ready  for  debugging: 

a.  Executive  Program 

b.  Interrupt  Process  Routines 

c.  Software/System  Development  (Debug  Aids)  Program 
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d.  Paper  Tape  Load/Dump  Routines 

e.  Terminal  I/O  Routines 

f.  Initialization/Restart  Program. 

Development  plans  are  being  formulated  for  the  radio  interface  software  and  digital  hard¬ 
ware  test  software  elements. 


Section  2 


Pseudonoise  J^iread  Spectrum  Analysis 


2.  1  INTRODUCTION 


The  ARPA  experimental  repeater  utilizes  a  pseudonoisc  spread  spectrum  (PNSS)  signaling 
technique.  This  section  of  the  quarterly  report  is  a  tutorial  discussion  of  this  spread 
spectrum  approach.  In  the  experimental  repeater  minimum  shift  keying  (RISK),  which  can 
be  defined  as  phase  continuous  FSK  with  a  modulation  index  equal  to  0.5,  is  used  to  modulate 
the  carrier.  The  reasons  for  this  choice  are  delineated  in  the  Experimental  Repeater 
Design  Plan  The  more  classic  modulation  technique,  PSK,  is  used  throughout  this  section 
of  the  report.  This  was  done  to  enable  the  reader  to  concentrate  on  the  spectrum  spreading 
discussion  itself  and  not  be  encumbered  with  the  additional  conceptual  complexities  of  RISK. 
The  fundamental  characteristics  and  properties  of  PNSS  signaling  are  essentially  the  same 
regardless  of  the  exact  form  of  angle  modulating  the  carrier. 


Some  of  the  principal  advantages  of  the  spread  spectrum  approach  discussed  and  spread 
spectrum  techniques  in  general  are  their  operation  in  the  presence  of  strong  interfering 
signals  such  as  deliberate  jamming  and  other  users.  In  addition,  the  signals  are  not  easily 
det jctable  because  they  tend  to  blend  into  the  background  noise.  Another  advantage  is  that 
the  spread  spectrum  signals  can  aid  in  multipath  rejection.  It  is  relatively  easy  to  also 
Kicorporate  addressing  into  a  spread  spectrum  system.  Historically,  one  principal  dis- 
advantage  has  been  the  period  of  time  required  for  synchronization.  Excessive  synchronization 
time  is  particularly  inconvenient  in  a  dynamic  interactive  type  of  environment.  Rapid 
advances  in  technology  have  recently  uncovered  some  exciting  solutions  to  this  problem. 


This  section  of  the  report  begins  with  a  basic  discussion  of  PNSS.  The  generation  and  basic 
properties  of  pseudorandom  sequences  are  then  discussed  in  the  context  of  PNSS  systems. 
Pseudonoise  code  synchronization  techniques  are  considered  and  a  quantitative  analysis  is 
given  developing  the  processing  gain  and  multiple  access  capabilities  of  PN  spectrum 
spreading  over  a  linear  channel. 


2.2  PSEUDONOISE  (DIRECT  SEQUENCE)  SPREAD  SPECTRUM 


Pseudonoise  (PN),  sometimes  called  direct  sequence,  is  a  well  known  and  widely  used 
spread  spectrum  signaling  method.  The  basic  form  of  a  PN  signal  is  that  produced  by  a 
pseudorandom ly  PSK  modulated  carrier.  In  this  modulation,  each  bit  of  information  is 
data  modulated  (eg,  +  1  in  binary  modulation)  upon  a  sequence  of  rf  carrier  chips  whose 
phases  are  modulated  according  to  the  output  of  a  PN  sequence  generator.  Ehe  PN  chip 
rate,  R  is  much  greater  than  the  data  rate.  A  generated  block  diagram  of  a  PN  spread 

spectrum  transmitter  is  shown  ir  figure  2-1.  A  typical  biphase  PN  signal  is  shown  in 
figure  2-2.  Although  PSK  is  mentioned,  the  characteristic  of  PN  spread  spectrum  can  be 
achieved  by  other  pseudorandomly  angle-modulated  rf  carriers.  Fundamental  to  the  PN 
modulation  is  the  concept  of  bandwidth  spreading  (the  transmission  of  signals  whose  rf  band¬ 
width  is  much  greater  than  the  message  bandwidth).  '"Tie  processing  gain  of  a  correlation 
receiver  is  equal  to  TW,  where  T  is  the  information  bit  time  and  W  is  the  spread  spectrum 
bandwidth. 
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Figure  2-1.  PN  Spread  Spectrum  Transmitter  Functional  Block  Diagram. 


Figures  2-3(a),  (b),  and  (c)  illustrate  the  spectrum  of  a  PN  modulated  signal.  The  original 
message  spectrum  is  a  narrow  spectrum  around  a  carrier  frequency,  with  a  main  lobe 
bandwidth  (null  to  null)  of  2R  for  PSK,  where  R  =  l/T  is  the  data  rate,  as  shown  in  fig¬ 
ure  2-3{b).  Since  the  PN  sequence  resembles  very  closely  a  pure  random  sequence,  the 
PN  modulated  signal  basically  occupies  the  same  bandwidth  as  that  of  one  rf  chip.  Hence, 
the  spectrum  of  the  PN  modulated  signal  has  a  main  lobe  bandwidth  (null  to  null)  of  R^  whereRc 
is  the  chip  rate  of  PN  modulation;  and  has  side  lobes  of  null-to-null  bandwidths 
equal  to  the  chip  rate,  as  shown  in  figure  2-3(c).  In  fact,  the  power  spectrum  of  the  PN 
signal  with  PN  code  period  K  at  chip  rate  Rc  has  an  envelope  of  the  form 

S/K  [sin  (jrf/R  )  /  (*f/R  )]2  for  PSK,  as  shown  in  figure  2-5,  where  S  is  the  carrier  power, 
c  c 

Thus,  the  total  available  rf  power  S  is  spread  over  a  wideband,  its  bandwidth  spread  is 
directly  proportional  to  the  chip  rate.  The  powe  r  level  at  each  harmonic  is  at  least  reduced 
by  the  factor  l/K  from  the  total  carrier  power,  for  PN  codes  with  period  K.  For  example, 
if  the  PN  chip  rate  is  10  MHz,  the  spectrum  will  have  a  main  lobe  with  20  MHz  bandwidth, 
and  side  lobes  with  10  MHz  bandwidths.  If  PN  codes  of  period  127  are  used,  the  power 
spectrum  density  will  be  at  least  21  dB  below  the  total  carrier  power  at  every  harmonic 
over  a  wideband.  Thus,  a  wideband,  low-power  density  signal  that  has  statistical  properties 
similar  to  noise  is  the  result.  Because  of  this,  PN  spread  spectrum  signals  will  have 
little  effect  on  narrowband  receivers  except  for  certain  near-far  relationships.  For  this 
system,  a  casual  listener  is  denied  access  to  the  baseband  information  since  it  can  be 
recovered  from  the  wideband  signal  only  through  correlation  with  a  stored  reference 
sequence  in  the  receiver  which  is  an  exact  replica  of  the  original  encoding  sequence. 
Moreover,  substantial  processing  gain  against  narrowband  jamming  can  be  achieved  through 
band  spreading.  For  example,  at  data  rate  of  10  kHz  and  spread  bandwidth  of  10  MHz,  the 
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AJ  processing  gain  will  be  30  dB.  Foi  this  case,  the  narrowband  interference  can  be  20 
dB  above  the  transmitted  signal  level,  and  at  the  output  of  the  demodulator  one  obtains 
signal-to-noise  ratio  of  about  10  dB. 

Figure  2-4  shows  a  block  diagram  of  a  typical  PNSS  receiver.  ^seau'ence 

multiplied  bv  a  PN  modulated  carrier  that  has  the  same  code  as  the  incoming  PN  sequence. 
When^cquisUion  has  been  obtained,  the  PN  signal  carrier  is  aligned  in  time  to  within  one 
S  o?Zrefe«nce  PN  sequence.  The  if.  output  is  integrated  and  dumped.  In  other  words, 
a  correlation  is  performed  for  each  bit  of  data  with  the  know-i  PSK/PN  sequence.  The 
message  is  decoded  according  to  the  correlation.  A+l  will  he  decoded  if  the  correlation 

a  vo  nthprwisp  a  -1  will  be  decoded.  The  desired  signal  lines  up  with  the  reference 
PNcode-  thus  after  multiplying  with  its  own  rtplica  (or  negative  replica),  the  PN  modulation 
is  removed  leaving  the  baseband  information,  in  frequency  domain  language  the  product 
of  the  desired  signal  and  the  local  reference  has  a  collapsed  spectrum  ecluallnS  of  the 
original  messie Modulation.  The  product  of  the  reference  signal  (itself  wideband)  and 
narrowband  jammer  (CW  jammer)  will  have  a  wide  power  spectrum,  spreading  the  jammer 
Sower  over  a wide  band.  Smilarly.  the  products  of  the  reference  signal  and  other  wideband 
interferences  (including  noise  and  other  uncorrelated  PN  modulated  signals  due  to  other 
users  in  a  multiple  access  operation-the  latter  resembles  very  closely  wideband  noise) 
are  also  spread  over  a  vide  band  The ' 

Unt  TuTnariwband  HlteTof  no. Se  bandwidth  1/T.  At  the  output  of  the  integrate  and 

appears  *a^  the^noise  powe^-^gTving^ri^Ui'the'factoi^Ttt^as  thTprocessi^g^dn^ainst 
narrowband  jammer  and  uncorrelated  PN  signals.  Thus,  AJ  and  multiple  access  capabilities 

are  achieved. 

2.3  PSEUDORANDOM  SEQUENCES 

In  this  section  the  generation  and  basic  properties  of  pseudorandom  sequences  will  be 
discussed  See  figure  2-5  for  a  functional  block  diagram  of  the  10-MHz  Band  $>read  1 
Period  PN  Signal  Generation.  PN  sequences  are  periodic  binary  sequences  wh°se  statistica 
properties  resemble  pure  random  sequences.  Yet  they  are  deterministic,  le.  predictable. 

The  following  properties  define  a  periodic  binary  sequence  to  lie  PN: 

In  each  period  of  the  sequence  the  number  of  I's  differs  from  the  number  of  0's  by, 
at  most,  1. 

Amoig  the  runs  of  l's  and  0's  in  each  period,  one-half  the  runs  of  each  kind  are  of 
length  1  one-fourth  of  each  kind  are  of  length  2,  one-eighth  are  of  length  3, 
so  on  as  long  as  these  fractions  give  a  meaningful  number  of  runs. 

The  autocorrelation  function  of  the  binary-  sequence  must  be  as  illustrated  in  figure  2-6(b). 
in  Other  words  the  number  of  agreements  subtracted  from  the  number  of  disagreements 
Twh^n  a  perifd'of  fte Tequepoe  is  con, pared  to  a  cyclic  shift  of  itself)  should  equal  to  the 
period  K  if  the  number  of  shifts  is  zero,  and  equal  to  -1  if  otherwise. 

As  compared  to  a  pure  random  sequence,  whose  autocorrelation  function  is  shown  in  fig- 
arc  2  0(a)  fte  n,ain  difference  lies  in  the  periodic  nature  of  the  PN  sequence  as  opposed 


a. 


b. 


c. 
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Figure  2-4.  PN  Spread  Spectrum  Receiver  Functional  Block  Diagram. 


to  the  nonperiodic  nature  of  a  pure  random  sequence.  The  power  spectrum  of  a  PN 
modulated  signal,  given  by  the  Fourier  transform  of  the  normalized  autocorrelation  function, 
is  then  equal  to  the  discrete  spectrum, 


S(f)  =  «0  ♦  £ 

n=-oo 

n/o 


sin  irf/R 


c 


where  K  is  the  period  of  the  PN  sequence,  is  the  PN  chip  time,  and  R£  =  l/Tc  is  the 

chip  rate.  When  K  is  large,  this  spectrum  is  a  line  spectrum  at  multiples  of  l/KT  with 

2  C 

an  envelope  given  by  l/K  [sin  *rf/R  /jrf/R  ]  ,  as  shown  in  figure  2-7.  Thus,  the  spectrum 

c  c 

occupancy  of  a  PN  spread  spectrum  signal  is  exactly  known. 

It  has  been  known  for  some  time  that  maximum  length  linear  shift  sequences  are 
pseudorandom.  Generation  of  these  sequences  can  be  achieved  by  linear  feedback  shift 
registers,  a  block  diagram  of  which  is  shown  in  figure  2-5.  Binary  sequences  generated 
by  linear  feedback  shift  registers  depend  upon  the  feedback  connections.  An  n-stage  linear 
shift  register  may  be  represented  algebraically  by  a  polynomial  with  binary  coefficients 
that  are  added  and  multiplied  modulo  2.  The  coefficient  of  the  X*  term  of  the  polynomial 
representing  the  shift  register  is  1  or  0,  depending  on  whether  or  not  the  feedback  tap 
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»(X)  - 1  ♦  x2  +  x3  +  x5 


(a)  A  Five  Stage  Register  and 

its  Polynomial  Representation 


Tx  lo 


(b)  Functional  Block  Diagram  for  10-MHz  Band  Spread, 
127  Period  PN  Signal  Generation  Using 
Surface  Acoustic  Wave  Device  (SAWD) 


Figure  2-5.  10-MHz  Band  Spread,  127  Period  PN  Signal 
Generation  Functional  Block  Diagram. 
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|t|<  tc 

OTHERWISE 


Autocorrelation  Function  of 
a  Purely  Random  Sequence 


ONE  PERIOD  R(rl 


Autocorrelation  Function  of  a  PN  Sequence 


Autocorrelation  Function, 


Figure  2-6 
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to  the  kU’  stage  is  or  is  not  connected^  *  /^aU^uTnc”  tea  sequel  derated  by  linear 

5  rr-“-  " 1 “  r  ir-ss- — ■? 

shut  register  is  primitive;  ie,  it  divides  1  ♦  X2"-1  but  no  polynomial  of  the  form  1  +  X 

,  x,  ^  on  i  The  number  of  distinct  maximal  length  sequences  generated  by  an^~st^ 
shiff  re^ster  >%, where  *<K,  is  the  number  of  positive  integers  less  than  K  an 

relatively  prime  to  K.  For  example,  a  five-stage  register  vUU  have  *(2  -l)/5  =  6  pnm.  tve 
polynomials,  each  generating  a  maximal  sequence  of  period  2  -1  -  3  . 

Generation  of  a  PN/PSK  signal  can -Iso  be  ^B^ng a 

as  shown  in  figure  2-6(b).  In  L  is  CJ*se<  inDlAulses  at  bit  times  with  data  modulation 

for^uccessive  information  bits  as  they  arrive. 

2.4  PN  CODE  CROSSCORRELATION  AND  COEXISTENCE 

When  correlation  is  performed  between  the  received  signal  and  a  replica  reference  PN 
signal at The receive?,  the  output  will  be  a  sum  of  the  following  correlations. 

a.  Correlation  between  desired  signal  and  reference. 

b  Correlation  between  reference  signal  and  other  undesired  PN  signals  sent  by  co¬ 
existing  terminals  in  the  same  common  band. 

c.  Correlation  between  reference  signal  and  noise. 

d  Correlation  between  reference  signal  and  narrowband  interference  (jamming). 

When  synchronization  is  obtained,  a.  will  be  directly  proportional  to  the  energy -f  the 
desired  signal  in  one  bit  of  data.  For  zero  mean  noise  and  jamming  c.  and  d.  will  both 

b  thfthameiOT  ^fKS-ap  e  riod  o^the1?  code?  f  o  r  cmitiniM)  us  S<:odef  depending  o^n  whether 

SSsSfeassssSSSSSSS" 

such  orthoganality  between  different  user  pairs  is  not  possible. 
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More  critical  than  the  degradation  in  signal  detection  performance  is  synchronization.  The 
ideal  autocorrelation  property  of  a  PN  code  is  used  for  synchronization  (see  paragraph 
2.  5).  However,  when  crosscorrelation  between  different  PN  codes  is  large,  the  receiver 
bit  sync  loop  may  lock  onto  the  crosscor relation  peaks  between  local  reference  and  undesired 
PN  signals  rather  than  onto  the  autocorrelation  peaks  of  its  own  communication  link,  thus 
failing  to  acquire  and  maintain  synchronization.  Therefore,  the  successful  use  of  PNSS 
systems  in  multiplexing  applications  depends  upon  the  construction  of  large  families  of 
encoding  sequences  with  uniformly  low  crosscorrelation  values. 

One  obvious  way  to  construct  these  codes  is  by  exhaustive  search,  however,  a  deterministic 
construction  algorithm  has  also  been  discovered  by  Gold  ^  .  Gold's  approach  is  to  construct 
families  of  non-maximal  sequences  whose  crosscorrelation  functions  are  uniformly  small, 
and  whose  autocorrelation  functions  are  not  two-valued,  but  nevertheless,  the  out-of- 
phase  autocorrelation  functions  are  small  and  have  the  same  bound  as  the  crosscorrelation 
functions.  Gold  showed  that  there  are  "preferred  pairs'  of  maximal  sequences  of  period 

K  =  2n-l,  where  n  is  the  number  of  shift  register  stages,  whose  crosscorrelation  function 
is  uniformly  bounded  by: 


n  + 1 

K  |R(  t  )  |  <  2  2  +1 

for  odd  n 

n  +  2 

K  |R(  r  )  |  <  2  2  +1 

for  even  n 

For  large  n,  this  bound  is  approximately  equal  to: 


K|  R(t)|  <  2  =  2/T 

Gold  has  given  an  algorithm  to  construct  these  two  preferred  pairs  of  maximal  sequences. 
Let  f^X),  f2(X)  denote  the  primitive  polynomials  of  these  two  shift  registers.  Gold  showed 
that  a  linear  feedback  shift  register  corresponding  to  the  product  polynomial  f^(X)*  f2(X) 

will  generate  2n+l  different  (non-maximal)  sequences  each  of  period  2n-i  and  such  that 
the  crosscorrelation  function  of  any  pair  of  such  sequences  are  uniformly  bounded  by 


2<n+2/2>+l 
bound. 


Their  out-of -phase  autocorrelation  function  has  the  same  uniformly  small 


To  illustrate  this  construction  we  consider  the  case  n  =  7.  In  other  words,  we  want  to 

construct  a  family  of  sequences  of  periods  2  -  1  =  127,  whose  pairwise  crosscorrelation 

is  uniformly  bounded  by  2^n+*^2)+  I  =  17  in  magnitude  (normalized  value  17/127  ~  0. 134). 

This  is  compared  to  the  maximum  autocorrelation  valt  e  of  127  (normalized  value  1). 

For  a  7-stage  shift  register,  there  are  <f>(27-l)/7  =*  126/7  =«  18  distinct  maximal  sequences.  It  is 
known  that  the  mutual  crosscorrelation  function  among  pairs  can  be  as  high  as  41  (see 
figure  2-8  for  two  maximal  sequences  in  18  that  have  this  high  crosscorrelation).  Using 
Gold's  algorithm,  one  can  pick  two  preferred  maximal  sequences  among  the  18  whose 
crosscorrelation  function  is  bounded  by  17.  They  have  the  following  primitive  polynomials: 

fx(X)  =  1  +  X  +  X2+X3  +  X7 


I 


f2(X)  =  1  +  X  +  X2  +  X3  +  X4  +  X5  +  X7 
References  are  provided  at  the  end  of  the  section. 
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Then  the  product  polynomial,  corresponding  to  a  14-stage  siiift  register, 

fx(X)  •  f2(X)  =  1  +  X2  +  X6  +  X8  +  x11  +  x12  +  x14 

according  to  Gold's  proof,  will  generate  129  different  non-ma>dmal  sequences  of  period 
127.  The  crosscorrelation  function  of  any  pair  of  such  sequences  will  be  bounded  by  17. 

It  is  of  interest  to  compare  Gold's  upper  bound  to  the  crosscorrelation  of  purely  random 
binary  sequences.  For  sequences  of  length  127,  the  variance  of  their  crosscorrelation 

is  Vl27  =  11.2  (the  mean  is  zero).  Thus,  Gold's  code  has  crosscorrelation  with 
magnitudes  less  than  the  2<r  point,  but  greater  than  the  1  o  point  of  the  crosscorrelation 
function  between  purely  random  sequences. 

2.  5  SYNCHRONIZATION  CONSIDERATIONS 

Since  the  receiver  must  have  an  exactly  matching  replica  of  the  transmitted  PN  signal  to 
perform  correlation  detection,  the  receiver's  copy  of  the  PN  code  must  be  in  very  close 
time  synchronization  with  the  desired  transmitted  wave.  In  the  case  of  I  NSS,  this  means 
that  time  synchronization  must  be  achievec  to  within  fractions  of  a  PN  c!  ip.  For  example, 
if  a  10-MHz  spread  spectrum  bandwidth  is  used,  the  receiver  local  key  stream  generator 
must  be  locked  well  within  0. 1  micr  second  of  the  epoch  of  the  l’eceived  signals.  Thus,  the 
time  synchronization  requirement  presents  one  of  the  critical  aspects  of  PNSS. 

For  coherent  detection  of  the  signal,  frequency  and  phase  lock  are  required  in  addition  to 
PN  code  and  chip  timing  alignment.  Frequency  and  phase  lock  may  be  achieved  through  the 
use  of  Costas  phase-lock  loops,  or  the  phase  shift  may  be  removed  from  the  incomming 
signal  by  squaring  a  biphase  or  squaring  a  quadriphase  signal  twice.  The  resulting  multiplier 
signal  is  then  used  as  a  reference  in  establishing  a  coherent  carrier  for  demodulation. 

PN  code  alignment  and  chip  timing  synchronization  makes  use  of  the  ideal  autocorrelation 
function  of  the  PN  code,  mentioned  in  paragraph  2.3.  A  classical  method  of  PN  signal 
epoch  detection  is  to  sequentially  search  over  all  possible  clock  time  positions  until  a 
correct  position  is  located,  as  indicated  by  an  increase  in  the  measured  correlation  output 
in  magnitude. 

A  block  diagram  showing  this  search  implementation  is  shown  in  figure  2-9(a).  Noncoherent 
detection  is  assumed  in  this  implementation. 

If  the  correlation  can  be  performed  at  wave  propagation  speed,  then  the  instantaneous  value 
of  the  correlation  between  the  incoming  PN  signal  and  a  fixed  PN  replica  at  different 
positions  of  time/code  alignment  will  be  available.  Surface  acoustic  wave  devices  (SAWD) 
or  charge  transfer  devices  (CTD)  can  be  used  to  perform  this  function  as  shown  in  fig¬ 
ure  2-9(b).  Thus,  the  TN  code  epoch  can  be  detected  in  real  time,  avoiding  the  long 
synchronization  time  required  for  a  serial  search. 

For  a  multiple  user's  environment  however,  near-field  transmitters  may  have  much  higher 
power  than  the  desired  signal,  causing  the  out-of -phase  crosscorrelation  between  the 
reference  PN  signal  and  the  undesired  PN  signal  to  have  values  very  close  to  (or  even 
higher  than)  the  in-phase  correlation  value  between  the  desired  signal  and  the  reference. 

This  will  cause  serious  problems  in  PN  epoch  synchronization.  Gold  codes  mentioned  in 
paragraph  2.4  may  improve  the  situation  somewhat,  but  serious  performance  degradation 
may  still  exist  if  the  undesired  coexisting  signals  are  at  much  higher  levels  than  the  desired 
one. 
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(a)  Sequential  PN  Epoch  Detection  (Noncoherent) 
for  PN  Code  Alignment  and  Chip  Sync 


(b)  PN  Epoch  Detection  (Noncoherent)  in  Real-Time 
Using  Charge  Coupled  Device 


Figure  2-9.  PN  Signal  Epoeh  Detection  Block  Diagram. 
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Complementary  (Golay)  codes  do  provide  a  two-valued  autocorrelation  function  with  zero 
out-of-phase  autocorrelation.  Frank  codes  (polyphase  codes)  also  have  zero  valued 
out-of-phase  autocorrelation.  Crosscorrelation  between  different  pairs  of  Golay  codes  is 
discussed  in  section  3  of  this  report. 

2.6  PERFORMANCE  ANALYSIS  OVER  A  LINEAR  CHANNEL 

The  AJ  processing  gain  and  multiple  access  capabilities  of  PN  spectrum  spreading  have 
been  explained  intuitively  in  paragraph  2-2.  Some  quantitative  analysis  will  be  given  here. 
For  clarity  f  presentation,  most  mathematical  details  are  omitted.  Also,  we  shall  assume 
PN/PSK  and  me  channel  to  be  linear.  Channels  with  hard  limiting  will  be  discussed. 

Consider  that  a  signal  bit  of  duration  T  is  received.  The  received  signal  waveform  r(t) 
will  be  a  linear  sum  of  the  following: 


r(t)  =  s(t)  +  n(t)  +  J(t)  +  C(t) 
where 

s(t)  =  The  desired  PN  modulated  signal  with  power  S 

n(t)  =  Wideband  Gaussian  noise  with  one-sided  spectral  density  NQ 

J(t)  =  Interference  signal  (or  jammer)  with  power  J  and  effective  bandwidth  Bj 

C(t)  =  Clutter  due  to  other  transmitters  in  the  same  multiple  access  band;  these 

are  wide  bandwidth  processes  basically  uncorrelated  with  the  desired  signal. 


For  simplicity,  we  neglect  C(t)  foi  the  moment.  For  white  Gaussian  noise,  the  optimal 
receiver  will  be  the  matched  filter.  Thus,  the  receiver  performs  a  correlation  between 
the  received  signal  r(t)  and  the  reference  signal,  which  is  an  exact  replica  of  the  desired 
PN  signal  (except  for  the  data  modulation  ±1)  and  is  itself  wide  band.  Denoting  this 
reference  signal  by  S'(t),  the  correlation  X  performed  over  one  bit  T  will  be: 


X 


S'(t)  r(t)  dt 


S’(t>  s('0dt 


n(t)  S'(t)  dt  + 


J(t)  S'(t)  dt 


(1) 


where 

S'(t)  =  S(T-t). 
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Since  data  is  decoded  according  to  X  ,  the  statistics  X  (which  is  Gaussian)  determines 
the  probability  of  bit  error.  In  particular,  the  output  signal-to-noise  ratio,  SNR^,  whose 

square  root  is  the  normalized  mean  of  X  ,  ie, 

a  x2 
SNR  =  — — 

0  ’*2 


where  X  =  mean  of  X  ,  =  variance  of  X,  determines  the  probability  of  bit  error  through 

the  well  known  error  function. 

T 

The  mean  of  X  is  /  S’ (t)  S(t)  dt  since  both  noise  and  jammer  are  zero  mean.  The  variance 
is  a  sum  of  variances  of  the  noise  and  jammer  terms.  Thus, 


SNR  = 
0 


rT 

2 

J  S(t)  S'(t)  dt 

-O 

2 

'  rT 

E 

1  n(t)S'(t)  dt 

+  E 

/  J(t)  S»(t)  dt 

m/ 

O 

^  -J 

J 

L°  J 

(2) 


By  Schwartz's  inequality,  and  because  S'(t)  is  directly  proportional  S(t),  we  have: 


T 

/  S2(t)dt-  / 


S'  (t)  dt 


0 

r  rT 

2 

r  /*t  i 

2 

E 

I  n(t)  S'(t)  dt 

+  E 

/  J(t)  S'(t)  dt 

.o 

_o 

(3) 


The  power  spectrums  of  the  output  process  obtained  from  correlating  S' (t)  with  n(t) 

2  2 

and  J(t)  are  [ <f> (f)  |  NQ/2  and  |#(f)|  Pj(f),  respectively,  where  <j>(f)  is  the  spectrum  of 
the  reference  signal  (ie,  the  transfer  function  of  the  matched  filter),  NQ  is  the  one-sided 
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power  spectral  density  of  the  noise,  and  Pj(f)  is  the  power  spectral  density  of  the  interfering 

jammer.  Since  pow'er  of  a  signal  is  given  by  integrating  its  spectral  density,  and  by  using 
the  Parseral's  theorem,  equation  (3)  can  be  rewritten  as: 


SNR 

0 


S^t)  dt 


|0(f)|  2  Pj(f)  df 


|0(f)|2  df 


(4) 


2 

The  spectrum  0(f),  due  to  PN  modulation,  has  a  (sin  X/X)  envelope  with  first  nulls  at 

Rc  (figure  2-7)  and  effective  bandwidth  \\  =  2Rc<  In  general,  the  jammer  bandwidth 

Bt  is  <<W.  A  reduction  in  jammer  power  is  thus  seen  in  equation  (4).  Also,  because 
J  2 

of  the  (sin  X/X)  shaping  of  0(f),  the  worst  case  jammer  will  be  a  narrowband  jammer 
centered  around  carrier  frequency.  If  J  is  the  total  jammer  power  of  the  jammer,  then 
the  resulting  SNRQ  is: 

T 

s2(t)  dt 

SNRo  =  ^ -  (5) 

f  *  J/2Ilc 

The  factor  l/2  in  the  term  j/2Rp  is  due  to  the  phase  angle  difference  between  jammer 

and  reference,  which  is  slowly  varying  between  0  and  2*.  Notice  this  l/2  factor  may  be 
lost  in  the  case  of  biphase  PN  when  the  CW  jammer  is  phase-locked  with  the  reference.  In 
quadriphase  (or  polyphase  PN),  however,  this  factor  of  l/2  is  always  present,  which  says 
that  biphase  PN  could  have  a  3  dB  less  AJ  gain  then  quadriphase  PN  under  the  hypothetical 
case  that  phase  lock  was  achieved  between  jammer  and  the  reference. 


If  is  the  power  of  the  signal  to  be  received,  K  is  the  number  PN  chips  for  bit,  and 

R  is  the  chip  rite,  then  equation  (5)  can  be  rewritten  as: 
c 


SNR 

0 


2K  S1 


N  R  +  J 
0  c 


(6) 


showing  the  number  of  PN  chips  per  data  bit  is  AJ  processing  gain. 
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The  case  of  clutter  is  very  similar  to  wideband  noise.  If  there  are  M-l  other  active  users 

in  the  channel,  each  with  different  power  levels  S  m  =  2,  3,  •  •  • ,  M;  and  if  p  is  the 

^  st 

mean  square  (normalized)  crosscorrelation  between  the  m  user's  PN  code  and  the  1 
(the  desired  PN  signal),  it  can  be  shown  that  the  correlator  output  SNRq  will  be: 


2K  S, 


SNR  = 
0 


(7) 


NR  +  J  + 
0  c 


S  K  p 
m  m 


Using  Gold's  upper  bound  on  PN  code  crosscorrelations  (see  paragraph  2.4),  pm  will  be 
bounded  by: 


2 


4 

K 


and  equation  (7)  is  reduced  to 


2K  S 

SNR  > - 

0 

NR  +  J  + 

0  c 

K 

Equation  (8)  shows  that  the  AJ  processing  gain  is  V  .  Since  K  is  the  number  of  PN  chips 
per  data  bit,  K  is  equal  to  the  product  T  (bit  time)  x  W  (chip  rate)  as  earlier  mentioned. 

As  an  illustration  for  the  foregoing  analysis,  we  consider  the  following  examples: 

2.6.1  Example  1.  AJ  Capability  Analysis 

Suppose  jammer  power  is  60  dB  above  the  signal.  Neglecting  noise  and  multiple  access, 
this  means  10  log10  (SNRQ)  =  10  log  2K  -  60.  Requiring  SNRQ  to  be  13  dB  (corresponding 

to  error  rate  4  x  10-6),  K  must  be  at  least  1  x  107.  In  other  words,  a  PNSS  system 
requires  a  chip  rate  in  excess  of  5  mega  chips  per  data  bit  to  achieve  70  dB  AJ  processing 
gain.  Of  course,  this  assumes  that  neither  notch  filter  nor  any  other  nonlinear  processing 
is  used. 


I 

I 
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2.  6.  2  Example  2.  Multiple  Access  Analysis 

Suppose  127  PN  chips  are  used  per  data  bit.  Assume  that  each  terminal  transmits  the 
same  power.  Neglecting  noise  and  jamming,  SNR^  =*  2K/4  (M-l).  If  the  output  SNR^  is 

required  to  be  10  dB,  then  M  =  8.  In  other  words,  no  more  than  8  simultaneous  users 
should  be  used.  In  this  analysis  Gold's  upper  bound  is  used  in  evaluating  the  code  cross- 
corj  elation.  Thus,  this  conclusion  may  be  somewhat  conservative.  Exact  calculations  can 
be  performed  once  the  code  sets  are  defined. 

2.  7  EFFECT  OF  HARD  LIMITING 

The  effect  of  hard  limited  channels  will  be  briefly  discussed.  Whenever  two  or  more 
independent  signals  pass  through  a  hard  limiting  repeater,  the  power  ratio  of  any  two 
component  signals  at  the  output  will  generally  differ  from  the  power  ratio  at  the  input,  the 
change  favoring  the  stronger  signal  at  the  expense  of  the  weaker  one.  The  worst  case  of 
suppression  involves  two  constant  envelope  signals,  one  greater  than  the  other  by  a  factor 
of  4  or  more.  The  relative  output  power  of  the  smaller  signal  is  reduced  by  another  factor 
of  4  (6  dB)  below  its  relative  power  at  the  repeater  input.  CW  jamming  is  one  example  of 
this  case.  In  fact,  it  is  known  that,  for  one  quadriphase  PN  user  and  a  tone  jammer  of  power 
much  greater  than  the  desired  signal,  the  presence  of  a  limiter  degrades  the  AJ  processing 
gain  predicted  in  paragraph  2.G  by  6  dB,  independently  of  the  jammer  phase  angle.  However, 
for  biphase  PN,  the  jammer  phase  angle  does  affect  performance.  In  fact,  the  biphase  PN 
signal  can  be  significantly  suppressed  by  the  tone  jammer  if  the  jammer  is  in  phase  with  the 
carrier  of  the  desired  signal. 

In  contrast  to  this  worst  case  (6-dB  degradation),  if  the  larger  of  the  two  signals  can  be 
characterized  as  a  Gaussian  process,  then  the  smaller  signal  is  reduced  by,  at  most,  1 
dB  below  its  input  relative  power.  When  any  large  constant-envelope  signals  arrive 
simultaneously  at  the  repeater  input  with  nearly  equal  power,  the  effect  is  like  the  Gaussian 
signal  case,  so  that  any  one  constant  envelope  signal  among  them  would  be  suppressed  by 
about  I  dB.  In  other  words,  without  narrowband  jamming,  just  in  the  presence  of  Gaussian 
noise  and  clutter,  the  output  SNRq  will  be  degraded  by  about  1  dB  from  that  predicted  in 

paragraph  2.  6  if  hard  limiting  occurs  in  the  repeater. 

2.8  NEAR- FAR  PROBLEM 

For  a  radio  network,  the  number  of  simultaneous  users  may  be  large  and  some  of  the  users 
may  be  located  in  the  vicinity  of  the  receiver,  and  thus,  having  much  higher  power  than  the 
desired  signal.  Coexistence  is  thus  a  critical  aspect  of  the  PNSS  radio  network.  This  pro¬ 
blem  will  be  solved  if  coexisting  users  use  completely  orthogonal  PN  codes.  Orthogonality, 
however,  requires  that  all  users  be  synchronized  and  is  not  practical  for  most  radio  net¬ 
works.  As  paragraph  2.6.2  illustrates,  to  guarantee  an  SNR^  of  more  than  10  dB  for  PNSS 

with  127  PN  chips  per  data  bit,  no  more  than  eight  simultaneous  users  should  be  active  at 
the  same  time.  This  means  that  the  ratio  of  total  undesired  signal  power  to  desired  signal 
power  should  not  be  more  than  8.2  dB.  This  figure  may  be  conservative,  but  it  illustrates 
the  drawbacks  of  a  constant  envelope  PNSS  system  in  terms  of  near-field  coexistence. 

Hybrid  systems  like  frequency  hopping  PN,  time  hopping  PN  may  improve  performance  in 
multiple  users'  environme.it. 
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Requiring  each  user  to  transmit  signals  at  lowest  level  and  to  increase  power  only  when  the 
repeater  failed  to  acknowledge  receiving  of  the  signal  may  improve  performance  to  some 
extent.  In  such  a  situation,  the  SNR^  at  the  ou'put  of  receiver  correlator  is  predicted  by 

equation  (8). 

Another  alternative  to  combat  the  near-far  problem  for  a  pure  PNSS  system  is  to  send  more 
chips  per  bit  for  a  distant  (low  power)  user  and  less  chips  per  bit  for  closer  (higher  power) 
users.  An  example  is  given  to  illustrate  this  application. 

Example: 

Consider  the  case  M  =  4  users,  two  of  which  operate  with  high  power  and  two  of  which 
operate  with  low  power. 

sx  =  10 

s2  =  10 

53  =  1 

54  =  1 

Assume  noise  is  negligible  and  jammer  power  J  =*  50.  Suppose  SNR^  is  required  to  be  10 
for  each  of  the  four  users,  and  solve  for  the  minimum  allowable  K2>  Kg  and  K4>  which 

indicate  the  allowable  information  rates  of  the  four  users.  From  equation  (8),  the  results 
are 


Kl  =  K2=49 

K  =  K  *  670 
3  4 

Since  the  chip  rates  are  constrained  to  be  equal,  the  information  rates  of  the  two  pair  of 
users  are  in  the  ratio  670/49  =  13.7.  That  is,  users  1  and  2  may  operate  at  a  rate  13.7 
times  as  great  as  users  3  and  4. 
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Code  Selection  for  the  Packet  Radio  Experimental  System 


3.  1  INTRODUCTION 

The  experimental  system  is  a  code  spread  system  where  each  kiti^ ‘  ~i'ded  int°  a  P^e“ 
determined  sequence  of  chips.  The  system  uses  MSK  (Minimum  Shift  Keying)  for  chip 
modulation.  The  MSK  chip  modulation  essentially  involves  transmitting  two  sets  of  ^phase 
modulate  1  chips,  staggered  by  a  chip  interval,  on  two  quadrature  subchannels  The  signals 
on  both  channels  are  cosine  w  eighted  over  every  2-chip  interval.  Therefore  the  resu  .  n 
transmitted  signal  is  a  phase  continuous  signal  despite  1  and  0  transitions  of  the  chip 
sequences  on  each  subchannel. 

Hypothetically  if  the  two  quadrature  subchannels  would  be  independent  channels  then  the 
chip  sequence  on  either  channel  could  conceivably  be  selected  from  any  suitable  f ami  y 
of  codes  to  give  the  best  autocorrelation  peaks  with  minimum  side  lobes.  Practically, 
however  the  two  subchannels,  X  and  Y,  are  not  isolated.  Consequently,  through  individual 
subchannel  matched  filters,  the  MSK  signal  when  correlated  at  the  receiver  results  in 
crosscorrelation  components  between  subchannels  in  addition  to  the  normal  autocorrelation 
components.  This  is  shown  by  figure  3-1.  Therefore,  the  codes  on  channels  X  and  Y  should 
have  minimum  crosscorrelation  side  lobes  in  addition  to  the  normal  requirements  of 
minimum  autocorrelation  side  lobes.  It  has  to  be  recognized,  however,  that  since  the 
signals  between  two  channels  are  staggered  by  a  chip,  the  baseband  crosscorrelation 
terms  do  not  simply  add  to  the  baseband  autocorrelation  terms.  In  fact,  the  cosine  weighting 
over  2-chip  interval  changes  this  even  further.  It  has  been  derived  in  Packet  Radio 
Temporary  Note  No.  93  that  the  magnitude  of  the  receiver  matched  filter  response  is 

given  by  the  following  equation: 


Envelope  Magnitude  I 


T  ) 
c 


RYY(t  " 


T  )] 
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2T  )  -  R 
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The  above  equation  suggests  that  the  baseband  correlation  terms  of  Rxx  and  Ryy  are  in 
phase,  but  Ryx  lags  them  by  a  chip,  and  R^  leads  them  by  a  chip.  A  computer  program, 

MSKCOR  for  the  MSK  baseband  correlation  is  developed  in  Packet  Radio  Temporary 
Note  No.  93.  The  results  of  the  model  described  tnerein,  have  been  confirmed  using 
EFT  (Fast  Fourier  Transform). 


It  is  apparent  from  the  above  matched  filter  respmise  equation  that  if  Ryx  -*  RYY  go  to 

zero  everywhere  except  for  the  perfect  match,  then  the  problem  reduces  to  one  of 
minimizing  the  difference  of  crosscorrelation  terms  staggered  by  2  chips.  1  he  Golay 
complementary  pairs  of  sequences  have  the  property  that  R^  +  Ryy  equals  zero  except 


for  zero  shift.  But,  the  crosstalk  between  the  two  subchannels  needs  to  be  investigated 
for  the  MSK  application. 


3-1 


code  selection  foT’  the  packet  radio  experimental  system 


TRANSMITTER 


RECEIVER 


Figure  3-1.  Crosscorrelation  Components  Between  Subchannels. 


The  packet  radio  experimental  system  uses  dual  data  rates  of  100  kb/s  and  400  kb/s.  This 
imposes  additional  constraints  on  any  code  selection  scheme.  The  codes  for  the  two  rates 
must  be  different,  and  the  crosscorrelation  between  codes  must  be  minimum  to  avoid 
interference  on  one  data  rate  detection  interfered  by  transmission  of  the  other  data  rate 
signal.  The  experimental  sjstem  will  have  a  spread  of  128  chips/bit  for  the  100-kb/s  rate 
and  a  spread  of  32  chips/bit  for  the  400-kb/s  rate,  thereby  resulting  into  12.8  megachips/’s 
for  both  data  rates.  There  are  5.25  if.  cycles  per  chip  to  result  in  the  first  stage  if.  fre¬ 
quency  of  67.27  MHz.  This  is  the  frequency  o'  the  transmit  and  the  receive  surface 
acoustic  wave  devices  (SAWD). 

The  following  paragraphs  present  the  solution  to  generating  the  desired  length  Golay  codes, 
and  discuss  the  code  selection  process  on  the  basis  of  minimizing  the  autocorrelation  and 
the  crosscorrelation  side  lobes.  It  must  be  emphasized  that  because  of  the  baseband  model 
approach  developed  in  Packet  Radio  Temporary  Note  No.  93  a  viable  code  optimization 
scheme  has  been  possible.  For  example,  due  to  the  referenced  approach,  5760  possible 
distinct  Golay  pairs  of  64  length  (required  for  the  low  data  rate  128-chip  spread  case)  are 
generated  and  autocorrelated  in  approximately  45  minutes  on  the  Univac  1108. 

Paragraph  3.2  describes  the  Golay  code  generation  algorithms.  Paragraph  3.3  discusses 
the  GOLSEL  computer  program  designed  to  generate  the  Golay  codes,  and  to  optimize 
their  selection  process  for  the  desired  autocorrelation  characteristics  for  single  data 
rate  systems.  Paragraph  3.4  provides  the  results  of  code  autocorrelation  of  the  16-length 
Golay  pairs,  and  tne  64-length  Golay  pairs  for  the  experimental  system  code  investigation. 
Paragraph  3.  5  describes  the  use  of  MSKCOR  for  crosscorrelating  of  dual  data  rate  codes, 
and  also  provides  the  results  of  the  dual  data  rate  codes  crosscorrelation  for  the  ex¬ 
perimental  system.  Paragraphs  3.6  and  3.7  summarize  the  final  code  pairs  selected 
along  with  full  details  of  their  side  lobe  characteristics. 

3.2  GOLAY  SEQUENCE  GENERATION 

The  basic  mechanics  of  generating  sequences  of  length  2N  involve  appropriately  combining 
sequences  of  length  N.  Begin  with  the  known  kernel  of  length  2  and  build  up  longer  sequences 
with  lengths  of  4,  8,  16,  23,  64,  etc.  One  conjecture  is  that  for  a  sequence  length  of 

2  ,  there  are  n.' 2  (n  >  2)  possible  distinct  Golay  pairs^.  Based  on  this  result,  for  a 


Reference  is  found  at  the  end  of  the  section. 
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sequence  length  of  64,  n  =  6,  and  there  are  5760  distinct  pairs.  Similarly,  for  a  length  of 
16  (applicable  for  the  high  data  rate  32-chip  spread  case)  there  are  48  distinct  pairs  of 
sequences. 

If  X  and  Y  denote  complementary  sequences  subscripted  by  the  desired  length,  then  a 
particular  set  of  algorithms  can  be  used  to  generate  the  predicted  number  of  distinct 
sequences  and  their  Golay  complements.  The  algorithms  are  presented  in  table  3-1. 

3.3  CODE  SELECTION  PROGRAM 

A  GOLSEL  program  is  written  to  select  a  given  length  code  based  on  certain  predefined 
constraints.  The  program  includes  generation  of  codes  based  on  table  3-1.  Figure  3-2 
is  a  representative  flow  chart  of  the  mechanics  of  the  optimization  process. 

Basically  as  is  apparent  from  figure  3-2  also,  the  desired  code  length  is  called  as  an  input 
to  the  program.  The  length  has  to  be  2n  (where  n  >  2)  for  this  particular  program  because 
the  generation  is  based  on  a  kernel  length  of  2.  The  starting  code  is  a  pair  of  length  4 
where  the  code  X  is  1  1  1  0  and  the  code  Y4  is  1  1  0  1.  Higher  length  codes  are  then 

generated  through  subroutine  GOLGEN  by  means  of  algorithms  of  table  3-1.  First,  the 
length  8  pairs  are  generated.  As  pointed  out  in  note  5  of  the  referenced  table,  the  distinct 
six  pairs  of  length  8  are  generated  by  using  algorithms  #1-6.  These  algorithms  are 
referred  to  as  For  example,  if  M^  =  5,  then 

Xg  =  x4  y4  =  (1  1  1  0)1(1  10  1)  =  11111001 


Y  =  X  Y  =  (1  1  1  0)  (0  010)  =  10101100 

o  ^  1  ^  ^ 

That  is,  the  code  pairs  X4  and  Y4>  or  X4  and  Y4  are  simply  interleaved. 

Further,  to  obtain  the  48  distinct  pairs  of  length  16  for  each  M^,  algorithms  #1-8  must 
be  used.  These  algorithms  are  referred  to  as  Mg.  For  example,  if  Mx  =  5  and  Mg  =  3,  then 

X  =X*  !  Y  =  (1111100  1)*  i  (10  10  110  0) 

lo  o  o 

=  1001111110101100 
Y.c  =  x*  |  Y„  =(1111100  1)*  i  (10101100) 

lo  8  '  o 

=  1001111101010011 

Although  for  convenience,  the  flow  chart  in  figure  3-1  applies  to  generating  and  optimizing 
codes  up  to  length  16,  the  generation  and  optimization  of  higher  length  codes  are  an 
extension  of  this  procedure.  It  may  be  noted  for  further  reading  that  Mg  will  refer  to 

algorithm  assignments  for  extending  length  16  to  length  32,  and  M4  will  refer  to  algorithm 

assignments  relating  to  extension  of  length  from  32  to  64.  M^  can  vary  from  1  to  10,  and 

M  from  1  to  12.  Therefore,  the  maximum  number  of  distinct  64  length  Golay  pairs  is 
4 

M.  M  M  M  =  6.  8.  10. 12  =  5760. 

X  ^  o  ■ i 
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Table  3-1.  Algorithms. 


2W  N\2^  N 

•  X2N  *  V2/XN 

•  X2N  =  V2/YN 

°*  X2N_YN^4/XN 
L  X2N  =  V8/YN 
2*  X2N  =  V8/XN 


2N  N^2/ 
Y2N  =  YN 


Y2N  V 


Y  =  Y 

2N  n^4/ 

Y  =  v 
2N“ 

Y  =  Y 

2N  N^8/ 


XN  °r  YN  imply  secluence  inverted;  i.  e. ,  if  X  is  1  1  1  0,  then 
means  0  0  0  1. 

*  * 

XN  or  YN  imply  se<luence  conjugate;  i.  e. ,  if  X4  is  1  1  1  0,  then 
* 

X4  means  0  111. 

XN  !  Yn  imply  a  union  two  sequences.  If  X4  is  1  1  1  1  and 
Y .  is  0  0  0  0,  then  XXT  !  YXT  means  1111000  0. 
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Notes: 


4. 


Table  3-1.  Algorithms  (Cont). 


Y  implies  interleaving  two  sequences  every  I;  i.  e. ,  if  X  is  1  1  1  1 

✓  N 


2^4 


N^l/ 

and  Y  is  0  0  0  0,  then  X  Y  means  10  10  10  10,  and  X4 

4  4 

means  1  1  0  0  1  1  0  0  and  so  on. 

To  obtain  X4  and  Y4>  use  algorithm  1.  To  obtain  the  six  distinct  X8,g 

and  Y  ,  ,  use  algorithms  1  through  6.  To  obtain  the  48distinct  X^G,g 

and  Y  use  algorithms  1  through  8  to  each  of  the  distinct  Xg(g  and 

Y  ,  .  To  obtain  the  480  distinct  X^.g  and  Y32,s,  use  algorithms  1 

through  10  on  each  of  the  distinct  X1G,g  and  Y1(.,g.  To  obtain  the  o7o0 

distinct  X,,.,  and  Y,..lo,  use  algorithms  1  through  12  on  each  distinct 
04'S  s 

X32  andY32,  etc. 


It  is  evident  from  figure  3-2  that  as  each  code  of  a  desired  length  is  generated  ^  first  goes 
through  the  MSKCOR  subroutine  before  the  next  distinct  code  is  generate  . 
subroutine  uses  the  MSK  correlation  model  approach  described  in  Packet  Radio  Temporary 
Note  No  93  The  correlation  process  normally  involves  correlating  2  bits  propagating  past 
fbu  that  is  either  CODE  followed  by  CODE  can  be  correlatedwith  the  same  CODE  to  give 
autocorrelation  in  the  sequential  mode,  or  CODE  followed  by  CODE  can  be  correlated  w.th 
‘the  CODE.  In  later  presentations,  the  first  mode  of  correlation  will  be  normally  ^red 
to  as  C  C  *  C  and  the  second  mode  of  correlation  will  be  referred  to  as  C  C  C. 

MSKCOR  subroutine  is  preselected  for  either  of  the  two  modes  of  correlation,  that  is, 

CC*CorCC*C. 

The  code  first  propagating  past  the  matched  filter  (or  the  first  bit  propagating  past  the 
detector  SAWDS)  results  in  autocorrelation  side  lobes  which  are  ignored.  inn  y, 
side  lobes  resulting  at  the  tail  end,  when  the  code  is  not  followed  by  anything,  are  als 
torcd  Th^  onS  relevant  side  lobes  considered  are  the  ones  that  reflect  a  norma  mode 
of  data  detection  where  <  ach  data  bit,  except  for  the  first  and  the  last  is  erttierjo  owe 
or  precede  1  by  another  oata  bit.  The  mean,  variance,  and  the  peak  of  the  relevant  side 
lobes  are  written  and  sto  -ed  for  each  distinct  pair  of  Mj.  and  Mg  for  the  10-length  Golay 

•  Tho  mpnn  Deak  and  variance  of  the  next  code  can  be  compared  for  any  desired 
PoplfmiZa^on  sSegPy  selected.  H»e  present  GOLSEL  configuration  is  based  on  listing  mean 
variance  and  peak  of  all  codes,  and  finding  the  code  that  gives  the  minimum  peak  and  the 
minimum  mean.  Equal  weighting  is  given  to  all  of  the  parameters.  In  fact  die  aPProac 
is  presented  here  to  demonstrate  that  the  optimization  strategy  can  be  predetermined  and 
implemented  accordingly;  however,  the  result  would  be  more  important  for  single  data 
rate  detector  cases  since  in  the  case  of  the  multiple  data  rate  detectors  the  crosscorrelation 
side  lobes  also  must  be  considered.  A  typical  output  from  the  GOLSEL  program  is  listed 
in  table  3-2  for  reference  for  the  16-length  Golay  pairs. 

Note  that  in  order  to  minimize  the  computer  processing  time  all  of  J;he  j7^°  c°des  of 
length  64  need  not  be  subjected  to  GOLSEL  in  either  the  C  C  *  C  or  the  C  C  C  mode. 
Instead  only  the  C  *  C  correlation  needs  to  be  done  for  all  the  5760  codes.  A  gener 
observation  from  the  C  *  C,  the  C  C  *  C,  and  the  C  C  *  C  correlations  of  16-length  pairs 
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•K"  IS  THE  DESIREO  GOLAY  LENGTH 


K  ■  16 
P  0 
□  -  0 

MEAN  (11  -  100 
VAR  (II  100 
PEAK  (II  100 
HOOP  1 


OAT  A 

CX4  1110 
CY4  1101 


00  1000 
Ml  1  6 


CALL 

MSKCOR 
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Figure  3-2.  GOLSEL  Program  Flow  Chart. 
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is  that  the  minimum  peak  side  lobes  of  C  *  C  result  into  small  peak  side  lobes  for  C  C  *  C 
and  C  C  *  C  Then,  the  C  C  *  C  and  the  C  C  *  C  cai  -e  performed  on  a  much  smaller 
number  of  distinct  pairs  that  correspond  to  minimun  peak  side  lobes  for  C  *  C  This 
time  saving  approach  is  utilized  for  the  low  data  ratt  '128-chip  length  or  64-longth  pairs) 
codes  for  the  experimental  repeater. 

3.  4  AUTOCORRELATION  OF  CODES 

The  C  C  *  C  and  the  C  C  *  C  autocorrelation  results  of  all  of  the  16-length  and  the  64- 
length  codes  are  summarized  in  this  section.  The  16-length  code  pairs  are  applicable  for  the 
high  data  rate,  and  the  64-length  code  pairs  are  applicable  for  the  low  data  rate. 

Table  3-3  lists  the  C  C  *  C  and  the  C  C  *  C  autocorrelation  results  of  the  16-length  high 
data  rate  code  pairs.  As  discussed  in  the  earlier  paragraph  on  code  generation  algorithms, 
there  are  a  total  of  48  distinct  Golay  code  pairs  in  this  length  category.  varies  from  1 

through  6.  For  each  value  of  Mg  varies  from  1  through  8.  Henceforth,  the  distinct 

pairs  will  be  referred  to  as  (M  M2).  Assuming  that  a  greater  weighting  is  given  to 

minimizing  the  peak  of  all  side  lobes,  it  is  evident  from  J_able  3-3  that  the  following  13 
code  pairs  give  the  lowest  peaks  for  both  C  C  *  C  and  C  C  *  C: 

(1,  1),  (1,  2),  (1,  6),  (2,  8),  (3,  5),  (3,  7),  (5,  2),  (5,  3),  (5,  4),  (5,  7),  (5,  8), 

(6,  2),  (6,  7). 

The  above  observation  corresponds  to  the  minimum  peak  side  lobe  value  of  8  for  C  C  *  C 
and  6  for  C  C  *  C.  For  a  system  with  a  single  data  rate,  code  selection  with  minimum 
mean-variance  combination  could  be  taken  from  the  above  13  codes.  But,  since  the  ex¬ 
perimental  repeater  uses  dual  data  rates,  this  decision  is  subject  to  the  crosscorrelation 
side  lobes  between  the  two  data  rate  codes. 

Table  3-4  lists  the  C  C  *  C  and  the  C  C  *  C  autocorrelation  results  for  the  64-length  low 
data  rate  codes.  As  discussed  in  the  last  section,  not  all  of  the  5760  codes  are  subjected 
to  the  C  C  *  C  and  the  C  C  *  C  mode  of  autocorrelation.  Only  41  out  of  the  5760  codes  that 
result  into  a  minimum  peak  of  13  in  a  nonrepetitive  code  autocorre’ation  mode  (simply 
C  *  C)  are  considered.  The  C  C  *  C  and  the  C  C  *  C  autocorrelation  results  for  these  41 
codes  are  summarized  in  table  3-3.  M  Mg,  Mg,  M4  vary  from  l  to  a  maximum  value 

of  6,  8,  10,  12  respectively.  The  distinct  pairs  of  this  64-length  class  are  later  referred 

to  as  (M  ,  M  M  ,  M  ).  There  are  three  codes  in  length  64  that  give  the  minimum  value 
12  3  4  _ 

of  peak  side  lobes  for  the  C  C  *  C  and  the  C  C  *  C  autocorrelation: 

(5,  7,  6,  11),  (6,  8,  6,  11),  (5,  3,  8,  3) 

As  presented  in  table  3-4,  for  each  of  the  above  three  Golay  code  pairs  of  64  length,  the 
peak  for  the  C  C  *  C  is  16,  and  the  peak  for  the  C  C  *  C  is  18.  There  are  codes  with  smaller 
peak  values  of  14  for  the  C  C  *  C  autocorrelation,  but  they  have  much  higher  peaks  for  the 
C  C  *  C  autocorrelation. 
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3.  5  CROSSCORRELATION  OF  DUAL  DATA  RATE  CODES 
The  13  codes  of  length  1G  and  the  3  codes of  length  84 .which 

based  on  the  minimum  peak  side  lobe  criterion,  ar  .  y  problem  as  opposed 

length  codes!  The  fol.owi’ng  paragraphs 

deal  with  these  problems. 

3.5.1  Crosscorrelation  of  High  Data  Rate  Codes  into  Low  Data  Rate  Codes 
a  low  data  rate  128-chip  spread  code;  that  is  the  inte“®re  there  is  4  bits  worth  of  high 

of  short  code  arriving  at  the  receiver  followed ^by  a  |£°«  C  °r  m  u  writte„  to  handle 
5-short  bit  length  s gnal  the  two  long  bit  length 

“;7r^uired,tSu!:Tas7taee  short  code  hfc.  of  the  signal  can  be  added  arbitrarily. 
Figure  3-3  shows  the  mechanics  of  this  scheme  of  crosscorrelation. 

a  u  •  r  0  o  o  thp  crosscorrelation  side  lobes  that  would  enter  consideration  are 

£:r;ss  — “  sss'SMrK 

—SrSs-SrSs 

S’S.'SEStSSSi This  i. ...  i»  ...i—  »•  «  “• 

high  data  rate  codes  into  the  low  data  rate  codes. 

the  high  data  rate  code  is  analyzed. 

3.  5. 2  Crosscorrelation  of  Low  Data  Rate  Codes  into  High  Data  Rate  Codes 

The  crosscorrelation  of  the  low  data  rate  long  codes  with  the  high  data  rate  short  codes  is 

refbVUeSof  128° chip6  C^e^efore^ithe ^su mp'tion^is  UiaT the Inte^e ring  code  signal 
Sfb? C  folfowed  * C ;  or  ctuowed  by  C.  The  process  need  not  be  repeated  for  C 
followed  by  C,  or  C  followed  by  C,  as  the  overall  results  will  be  identical. 


3-13 


code  selection  for  the  packet  radio  experimental _ system 


1  BIT  OF 

HIGH  DATA  RATE 
CODE 

hV  , 

I 


CCCCjCCCC 

cccccccc 

cccccccc 

cccccccc 

CCCC|CCCC 

cccccccc 
cccccccc 
cccccccc 
_c  ccc.c  c  c  c 
C  ~C  c|c  I  REPEAT  THE 
I  'ABOVE  16 
I  j  COMBINATIONS 


FOR  EACH  OF  THE  32 
COMBINATIONS  LISTED 


TO  THE  LEFT,  CONSIDER  ALL 
SIDE  LOBES  APPEARING 
WITHIN  THIS  INTERVAL 
IN  LENGTH  EQUAL  TO  A 
SHORT  CODE  BIT 


Figure  3-3.  High  Data  Rate  into  Low  Data  Rate  Crosscorrelation. 
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le  3-6.  Length  32  Code  (5,  3)  into  Length  128  Code  (6,  8,  6,  11)  Crosscorrelation  (Cont) 
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The  MSKCOR  program  is  set  up  to  correlate  a  2-bit  length  signal  with  a  1-bit  length  cor¬ 
relator  (matched  SAWD)  with  both  the  signal  and  the  correlator  bits  having  identical  number 
of  chips.  Therefore,  the  first  64  chips  of  the  incoming  signal  are  correlated.  As  the  first 
two  chips  exit  the  correlator,  another  two  chips  are  added  from  the  remaining  longer 
string  of  the  long  signal  and  again  subjected  to  the  MSKCOR  program.  Note  that  whenever 
64  chips  of  the  input  signal  (equivalent  to  twice  the  length  of  the  correlator)  are  correlated 
with  a  32-chip  coded  SAWD,  the  program  would  normally  print  out  97  chip  correlation  points 
of  the  envelope  magnitude  of  the  baseband  correlation  curve.  Note  that  one  extra  correlation 
point  occurs  because  one  chip  delay  on  one  of  the  correlator  quadrature  channels  takes  care 
of  one  chip  staggering  between  the  two  quadrature  channels  of  the  MSK  signal.  The  32  end 
points  out  of  a  total  of  97  correlation  points  are  ignored  as  they  correspond  to  the  condition 
of  the  64 -chip  signal  followed  by  no  signal. 

Therefore,  only  the  first  65  (97  minus  32)  correlation  points  are  printed  out.  As  the  two  extra 
chips  are  added  to  compensate  for  the  two  exiting  chips  (implying  that  a  new  pair  from  the  64- 
length  code  pair  is  added  for  one  pair  leaving),  the  two  new  correlation  points  are  printed 
for  each  new  arrival  of  two  chips.  There  are  188  distinct  correlation  points  in  addition  to 
the  65  already  described  which  are  pointed  out.  The  mean,  variance,  and  peak  are  com¬ 
pared  for  each  set  of  codes  being  correlated. 

Table  3-7  summarizes  the  crosscorrelation  of  the  3  long  codes  into  the  13  short  codes  under 
investigation.  Table  3-7  provides  a  typical  computer  listing  of  the  crosscorrelation  side 
lobes  because  of  a  preferred  set.  The  mean.  '  ariance,  and  the  peak  of  the  set  correlated 
are  printed  out  at  the  bottom  line  of  table  3-8. 

3.6  CODESELECTION 

Referencing  tables  3-5  and  3-7,  it  is  apparent  that  the  short  code  (5,  3)  and  the  long  code 
(6,  8,  6,  11)  result  in  low  crosscorrelation  side  lobes.  Although  there  are  other  codes 
under  each  crosscorrelation  category  that  give  slightly  lower  peak  values  of  side  lobes, 
these  codes  are  not  suitable  when  both  categories  of  crosscorrelation  are  combined;  that  is, 
when  crosscorrelation  of  the  low  data  rate  into  high  data  rate  and  the  high  data  rate  into  the 
low  data  rate  are  considered  simultaneously.  Therefore,  the  best  Golay  code  pairs  for 
the  MSK  application  in  a  dual  data  rate  system  are  the  following: 

High  Data  Rate  Code;  (5,  3) 

Low  Data  Rate  Code:  (6,  8,  6,  11) 

Table  3-9  lists,  for  reference,  the  (5,  3)  codejilong  with  its  C  C  *  C  mode  of  autocorrelation, 
and  table  3-10  lists  the  (5,  3)  code  with  its  C  C  *  C  mode  of  autocorrelation.  A  con¬ 
firmatory  plot  of  the  (5,  3)  code  autocorrelation  in  the  C  C  *  C  mode  using  FFT  is  pro¬ 
vided  in  figure  3-4  to  validate  the  baseband  results  using  the  model  baseband  approach  of 
Packet  Radio  Temporary  Note  No.  93.  Notice  that  the  envelope  of  this  plot  matches  with 
the  peak-to-peak  envelope  at  every  chip  as  predicted  by  the  referenced  model. 

Tables  3-11  and  3-12  list  the  C  C  *  C  and  the  C  C  *  C  mode  of  autocorrelation  for  the 
long  (6,  8,  6,  11)  code  selected. 


e  3-7.  Low  Data  Rate  Code  (Length  128)  into  High  Data  Rate  Code  (Length  32)  Crosscorrelation. 


Table  3-8.  Crosscorrelation,  Length  128  Code  (6,  8,  6,  11)  into  Length  32  Code  (5,  3)  (Cont). 
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Table  3-8.  Crosscorrelation,  Length  128  Code  (6,  8,  6,  11)  into  Length  32  Code  (5,  3)  (Cont). 


code  selection  for  the  packet  radio  experimental  system 


! 


I 


■o 


rv  «  «  C 


K  00-^*^00  —  O'  N  N  mm 


r 


i 


o-wn<*/\oNc^o-Nn  <• 


J 


3-24 


Table  3-S.  High  Data  Rate  (Length  32)  C  C  *  C  Correlation. 
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li.ble  3-10.  High  Data  Rate  (Length  32)  C  C  *  C  Correlation. 
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3.  7  SUMMARY  OF  RESULTS 

Based  on  the  data  derived  in  the  preceding  paragraphs,  the  following  important  results  are 
summarized: 

a.  High  Data  Rate  (32-Chip  Length) 

Peak  Autocorrelation  Side  lobe  (C  C  *  C)  =  8.00 

Mean  Autocorrelation  Side  lobe  (C  C  *  C)  =  2.63 

Variance  of  Autocorrelation  Side  lobe  (C  C  *  C)  =  5. 10 

Peak  Autocorrelation  Side  lobe  (C  C  *  C)  =  6.00 

Mean  Autocorrelation  Side  lobe  (C  C  *  C)~=  2.44 

Variance  of  Autocorrelation  Side  lobe  (C  C  *  C)  =4.35 

Peak  Crosscorrelation  Side  lobe 
(due  to  low  data  rate  code)  =  12.  59 

Mean  Crosscorrelation  Side  lobe 
(due  to  low  data  rate  code)  =  5. 14 

Variance  of  Crosscorrelation  Side  lobe 
(due  to  low  data  rate  code)  =8.82 

Ratio  of  Autocorrelation  Peak  to  the  Peak  Autocorrelation 
Side  lobe (C  C  *  C)  =  12.04  dB 

Ratio  of  Autocorrelation  Peak  to  the  Peak  Autocorrelation 
Side  lobe  (C  C  *  C)  =  14.  54  dB 

Ratio  of  Autocorrelation  Peak  to  the  Mean  Autocorrelation 
(C  C  *  C)  Side  lobes  =  21.70  dB 

Ratio  of  Autocorrelation  Peak  to  the  Mean  Autocorrelation 
(C  C  *  C)  Side  lobe  =  22.  36  dB 

Ratio  of  Autocorrelation  Peak  to  the  Peak  Crosscorrelation 
Side  lobe  =  8 .  1  dB 

Ratio  of  Autocorrelation  Peak  to  the  Mean  Crosscorrelation 
Side  lobe  =  15.  96  dB 

b.  Low  Data  Rate  (128-Chip  Length) 

Peak  Autocorrelation  Side  lobe  (C  C  *  C)  =  16.0 
Mean  Autocorrelation  Side  lobe  (C  C  *  C)  =  6.24 
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Variance  of  Autocorrelation  Side  lobe  (C  C  *  C)  =  24.09 

Peak  Autocorrelation  Side  lobe  (C  C  *  C)  =  18.0 

Mean  Autocorrelation  Side  lobe  (C  C  *  C)  =  5.  11 

Variance  of  Autocorrelation  Side  lobe  (C  C  *  C)  =  22.28 

Peak  Crosscorrelation  Side  lobe  (due  to  high  data  rate  code)  =  32.  9 

Mean  Crosscorrelation  Side  lobe  (oue  to  high  data  rate  code)  =9.5 

Variance  of  Crosscorrelation  Side  lobe  (due  to  high  data  rate  code)  =  56  56 

Ratio  of  Autocorrelation  Peak  to  the  Peak  Autocorrelation  Side  lobe 
(C  C  *  C)  =  18.06  dB 

Ratio  of  Autocorrelation  Peak  to  the  Peak  Autocorrelation  Side  lobe 
(C  C  *  C)  =  17.04  dB 

Ratio  of  Autocorrelation  Peak  to  the  Mean  Autocorrelation  Side  lobe 
(C  C  *  C)  =  26.24  dB 

Ratio  of  Autocorrelation  Peak  to  the  Mean  Autocorrelation  Side  lobe 
(C  C  *  C)  =  28.0  dB 

Ratio  of  Autocorrelation  Peak  to  the  Peak  Crosscorrelation  Side  lobe 
=  11.  8  dB 


Ratio  of  Autocorrelation  Peak  to  the  Mean  Crosscorrelation  Side 
lobe  =  22.6  dB 
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Section  4 


Capture  Model  Description 


4.1  INTRODUCTION 


Tne  term  "capture"  has  been  defined  to  be  a  property  of  a  receiver  which  measures  its 
capability  to  distinguish  one  signal  from  an  ensemble  of  several  signals.  The  network  simu¬ 
lation  model  developed  by  Network  Analysis  Corporation(N.  A.  C. )  assumes  zero  capture- 
le,  it  two  or  more  packets  impinge  on  receiver,  none  are  recovered.  The  receiver  being 
developed  by  Collins  for  the  experimental  system  does  have  a  measure  of  capture  capability 
and  the  purpose  of  this  note  is  to  describe  this  property.  The  intent  is  for  the  capture 
phenomenon  vO  be  added  to  the  N.  A.  C.  model  with  the  objective  of  providing  additional 
reality  to  their  simulations. 


ro  begin  a  description  of  the  capture  model,  it  is  necessary  to  first  characterize  the  receiver 
Insofar  as  capture  is  concerned,  there  are  four  basic  characteristics  that  describe  the  state 
of  a  receiver.  These  are: 


a.  Whether  the  receiver  is  in  the  receive  enable  mode 

b.  If  in  the  receive  enable  mode,  whether  the  receiver  is  in  a  synchronization  period  or  a 
packet  reception  period 


c.  The  number  of  packets  available  to  be  captured  at  the  receiver 

d.  The  specific  characteristics  of  these  packets. 


The  first  characteristic  only  indicates  whether  the  radio  is  in  receive  mode.  Since  it  is  a 
half-duplex  operation,  the  radio  can  be  in  either  of  two  states:  transmit  or  receive.  The 
second  characteristic  determines  whether  a  receiver  is  receiving  a  packet  or  looking  for  a 
packet  preamble  for  synchronization.  Once  a  packet  has  been  synchronized  to,  then  the 

sy™hro™™tion  wiU  not  oc™r  again  until  this  packet  has  been  completely  received. 
The  third  and  fourth  characteristics  define  how  many  packets  are  available  for  reception  and 
their  specific  nature.  The  parameters  that  completely  describe  a  packet  for  the  capture 
model  are  the  following:  F 


a.  Time  of  arrival  of  the  packet  at  the  receiver 

b.  Power  level  of  the  packet 


c.  The  packet  length 

d.  The  data  rate  of  the  packet. 

These  four  parameters  are  the  only  packet  related  parameters  required  to  determine  the 
successfu.  capture  of  a  packet. 
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d.  The  data  rate  of  the  packet. 
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4.2  ASSUMPTIONS 

Before  continuing  the  description,  it  seems  appropriate  to  list  the  assumptions  that  simplify 
the  model  and,  more  importantly,  focus  the  attention  on  the  objective  of  modeling  capture. 

The  assumptions  are: 

a.  The  distances  between  all  network  elements  are  known. 

b.  The  receiver  has  two  detectors:  one  with  data  rate  100  kb/s  and  another  at  400  kb/s 
with  processing  gains  of  128  and  32,  respectively,  that  operate  simultaneously. 

c.  No  multipath  phenomenon  is  considered. 

d.  No  other  sources  of  interference  other  than  Gaussian  noise  and  packets  are  considered. 

e.  The  time  that  the  receiver  is  transferred  to  receive  enable  is  known. 

The  description  of  the  capture  model  is  divided  into  three  parts.  The  first  describes  the 
calculations  of  packet  parameters.  The  second  presents  the  operation  of  the  synchronization 
mode,  and  the  third  gives  the  operation  of  the  reception  mode.  Some  suggestions  lor  imple¬ 
mentation  of  the  model  into  the  N.  A.C.  simulation  program  are  given;  however,  the  total 
implementation  is  not  included. 

4.3  PACKET  CHARACTERISTICS 

It  has  been  assumed  that  the  packet  originations  and  their  time  of  origination  are  known.  In 
the  N.  A.C.  model,  this  is  obtained  from  an  event  generator.  The  two  significant  independent 
variables  are  thus  the  coordinates  of  the  transmitter  of  a  paeket  and  the  time  of  transmission. 

F ro m  these  two  variables,  <  ne  can  then  determine  the  time  of  arrival  of  this  packet  to  anv 
element  of  the  network  and  the  signal  level  of  the  packet  at  all  the  elements  in  the  network. 

The  derived  signal  level  at  each  element's  receiver  is  determined  from  a  propagation  atten¬ 
uation  calculation.  This  calculation  is  made  assuming  line  of  sight  for  repeater-repeater 
links.  For  repeater-terminal  links,  the  propagation  attenuation  is  based  on  Okumura,  et 

al  ,  propagation  data  which  is  described  in  detail  in  the  appendix. 

The  time  of  arrival  is  based  simply  on  free  space  delay  of  propagation.  The  data  rate  is 
implied  in  the  fact  that  it  is  known  which  element  transmitted  and  the  intended  receiver.  That 
is,  terminal-repeater  and  repeater-terminal  transfers  are  100  Kb/s,  while  repeater-repeater 
transfers  are  400  Kb/s.  The  packet  length  is  another  known  parameter  obtained  from  the 
event  generator;  and  for  simplicity,  it  might  be  assumed  that  all  packets  are  of  fixed  length. 

4.4  SYNCHRONIZATION  MODE 

With  the  description  of  the  derived  packet  characteristics  complete,  the  question  "How  does 
the  receiver  sift  through  these  available  packets  with  varying  characteristics  and  identify 
a  specific  packet?"  can  be  addressed.  The  answer  to  this  question  begins  with  the  synchroni¬ 
zation  process. 

It  has  been  pointed  out  that  packets  have  a  minimum  amount  of  information  associated  with 
them  that  uni  uely  describes  their  state.  This  can  lx  viewed  as  a  multiple  dimensioned 
vector  defining  a  packet  with  the  following  components:  To,  From,  Arrival  time.  Power  level. 


The 


reference  is  given  at  the  end  of  the  appendix. 
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Length .  I  he  synchronization  process  is  the  alignment  of  the  timing  of  the  receiver  and  a 
particular  packet  and  the  subsequent  discrimination  by  the  receiver  against  all  other  packets. 
Itis  based  on  the  operation  of  the  bit  synchronization  circuitry,  which  works  in  the  following  way. 

The  bit  synchronization  circuit  consists  of  a  phase-lock  loop  with  a  coarse  and  a  fine  tracking 
ability.  When  a  detector  in  the  receiver  is  enabled,  its  bit  synchronization  circuit  is  free- 
running  and  generating  three  windows  defined  as  Advance,  Retard,  and  Sync.  The  envelope- 
detected  correlation  peaks  coming  out  of  the  matched  filters  are  summed  together  and  a 
determination  is  made  as  to  which  window  has  the  largest  signal  residing  in  its  interval.  If 
a  signal  is  stronger  by  a  factor  of  3  dll  than  any  other  packet,  the  loop  will  slew  in  the 
advance  or  retard  direction  to  center  the  synchronization  window  over  this  larger  signal.  If 
this  signal  continues  to  be  larger  than  any  other  combination  of  signals  found  in  (he  Advance 
or  Retard  window  for  a  period  of  10  bit  intervals,  then  the  in-lock  condition  is  invoked.  This 
condition  narrows  the  Sync  window  down  from  (JT  >  to  2T  (T  is  period  of  chip)  and  enables 
the  preamble  detector.  c  c  e 

Also,  the  stronger  signal  must  not  coincide  in  its  correlation  peak  with  other  packet  corre¬ 
lations  within  2T  Therefore,  the  time  of  arrivals  of  this  signal  and  all  other  packet 

arrivals  must  be  divisible  by  the  bit  interval,  and  the  remainder  of  this  division  must  be 
greater  than  2T  to  ensure  proper  discrimination  of  the  packets.  Physically,  this  means 

that  the  time  of  arrival  of  the  correlation  peaks  from  the  desired  packet  must  be  separated  by 
at  least  2T^  from  the  arrival  times  of  any  correlation  peaks  from  other  packets.  If  this  is 

not  the  case,  then  the  end  of  preamble  will  be  assumed  to  go  undetected.  The  bit  synchroni¬ 
zation  circuit  will  continue  to  lock  up  to  this  signal  as  long  as  it  is  the  larger  signal. 

It  must  be  remembered  that  the  receive  enable  condition  could  have  been  obtained  at  a  point 
in  the  middle  of  the  strongest  signal  packet  in  which  its  preamble  has  already  been  transmitted. 
Therefore,  the  end  of  preamble  will  not  be  detected.  After  a  time-out  period  of  38  bit 
intervals,  the  synchronization  circuit  will  unlock  and  will  return  to  the  hunting  mode.  This 
means  that  the  successful  detection  of  end  of  preamble  (POP)  requires  that  the  strongest 
signal  must  reach  the  lock  condition  being  greater  than  any  other  packet  within  35T(j  or 

less  of  its  initial  arrival  time  to  ensure  that  the  preamble  is  detected.  Another  simplifying 
assumption  can  be  made  that  only  packets  at  one  data  rate  compete  for  synchronization. 

This  means  that  100-kb/s  packets  are  compared  for  synchronization  into  100-kb/s  data 
detector  but  not  400-kb/s  packets  and  vice  versa.  Figure  4-1  shows  a  flow  chart  of  synchro¬ 
nization. 


4 . 5  PACKET  REC E PT ION 


if  it  is  assumed  that  the  receiver  has  synchronized  to  one  packet  and  the  detection  decision 
is  made  within  a  2T  window,  then  an  expression  of  the  signal-to-noise  ratio  out  of  the  matched 

filter  can  be  obtained.  This  relation  reflects  the  time  discrimination.  The  general  equation 

is: 

2KS. 

SNR  =  - 


1 


M 


NR  +  X  KS  pt 
°  c  m=2  m  '■ 


m 
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where 


K 


NR 
o  c 


m 

¥  2  = 
1,  m 


M 


Power  level  of  synchronized  packet  signal 

Spread  factor  (number  of  chips  per  bit)  for  this  given  detector 
Thermal  noise  power  (constant) 

Power  level  of  m^1  packet 

Normalized  erosscorrelation  value  in  2T  window  of  PN  sequence  ana 

and  S  PN  sequence 
m 

Number  of  available  packets  at  the  receiver 


The  power  levels  of  the  ensemble  of  packets  impinging  on  the  receiver  are  known,  along  with 

the  value  of  K  which  is  128  for  the  100-kb/s  data  rate  and  32  for  the  400-kb/s  data  rate, 

dependim!  on  the  S.  data  rate.  The  noise  power  is  constant  and  equal  to  N  R  adjusted  by 
1  o  c 

the  noise  figure  of  the  receiver  (8  dli).  The  normalized  crosscorrelation  coefficient  P^  m 
is  dependent  on  the  codes  and  their  time  displacement.  For  example, 


1  if  t  <  2T 


PLI>>  =  63 


(128) 


2  if  '  ^2Tc 


where 


LL 


=  Normalized  crosscorrelation  of  low  data  rate  code  with  itself  and  represents 
the  displacement  of  correlation  peaks. 


The  value  of  PT  T  out  of  synchronization  was  obtained  through  computer  simulation  of  the 

L  L 

selected  codes.  Other  crosscorrelation  combinations  are: 


2 

1  if  | 

PHII  <*> 

12 

(32)2 

„  2.  , 

35 

1  LH  (r) 

(M 

CO 

II 

“  T  , 

146 

PIIL  (r) 

(128) 

if  r  >2T 


for  all  r 


for  all  t 


where 


- Z 

Pjj^  is  crosscorrelation  of  high  data  rate  code  with  itself 
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is  crosscorrelation  of  low  data  rate  code  into  high  data  rate  matched  filter 
is  crosscorrelation  of  high  data  rate  code  into  low  data  rate  matched  filter 


The  determination  of  r  ,  which  is  only  important  in  P^j  and  P^^  ,  is  obtained  from 
the  arrival  times  of  packet  m  and  the  synchronized  packet.  The  relation  is: 

|l'm  -  tJ 

t  =  fractional  part - ^ - 

b 


where 


T^  =  10  microseconds  for  low  data  rate 
2.5  microseconds  for  high  data  rate 

T  =  Arrival  time  of  packet  m 
m 

T^  =  Arrival  time  of  packet  1  (synchronized) 

_9 

It  should  be  remembered  that  2T^  =  15G  x  10  seconds. 

This  SNIt  will  vary  as  events  take  place  that  impact  the  receiver  and  as  packets  terminate; 
therefore,  the  signal-to-noise  ratio  should  be  recalculated  each  time  this  happens  and  a 
weighted  average  obtained  based  on  the  time  between  events. 

The  average  signal-to-noise  ratio  is  then  used  to  obtain  the  probability  of  a  bit  error  (Pe). 

The  probability  of  bit  error  as  a  function  of  signal-to-noise  ratio  for  differentially  coherent 
MSK  is: 

P  =  (1/2)0  -SNR/2 
e 

and  is  shown  graphically  in  figure  4  -2. 

With  a  value  of  P  ,  one  can  then  determine  the  probability  of  success  of  a  packet  of  length  n 
bits  from  the  equation: 

P  =  (1  -  P  )n 
n  c' 

The  P  provides  the  measure  needed  to  determine  if  a  packet  is  successful  received  error 
n 

free.  A  Monte  Carlo  method  can  lx*  used  to  determine  the  success  or  failure  of  a  packet 
reception.  This  could  be  the  selecting  of  a  random  number  Ixtueen  0  and  1  from  a  uniform 
distribution  v'ith  the  test  Ixing  that  if  this  number,  X,  is  greater  than  I5n,  then  the  packet 

has  one  or  more  errors  and  is  not  successfully  received.  If  X  is  less  than  P  .  then  the 

n 

packet  has  no  errors  and  is  successful. 
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4.6  IMPLEMENTATION  SUGGESTIONS 

The  capture  phenomenon  has  been  described  as  it  takes  place  in  the  Collins  design.  This, 
however,  very  likely  will  not  be  an  acceptable  algorithm  for  the  computer  simulation  network 
model  due  to  its  complexity.  Therefore,  reality  will  probably  need  to  be  tempered  with 
practicality  of  implementation. 

One  suggestion  is  to  partition  the  network  such  that  not  all  element  transmissions  impact  on 
all  other  network  elements.  The  criteria  of  selection  might  be  the  elimination  of  packets 
from  elements  whose  mean  signal  level  approaches  the  thermal  noise  level.  The  signal  level 
calculation  through  the  propagation  model  and  the  synchronization  model  are  basic  to  the 

capture  ability;  however,  the  packet  reception  model  might  be  simplified  by  assuming  no 

- 

packets  are  synchronous  with  the  synchronization  packet  and  the  (r)  terms  are 

relatively  very  small.  This  means  that  the  terms  in  the  denominator  of  the  SNR  equation 
associated  with  other  packets  might  be  neglected  and  only  one  signal-to-noise  ratio  calculated 
for  the  packet  transmission. 
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Appendix 

Okumura  Path  Loss  Approximation 


The  following  describes  a  computer  algorithm  which  has  been  generated  to  incorporate 
Okumura's  path  loss  attenuation  data  into  the  capture  model  for  the  N.  A.  C.  packet  radio 
simulation  program.  For  this  program,  the  average  effective  transmission  heights  for 
a  repeater  and  a  terminal  were  assumed  to  be  45  meters  (148  feet)  and  3  meters  (9.  8  feet), 
respectively. 

Figures  A-l  and  A-2  have  been  included  to  illustrate  the  process.  Figure  A-l  gives  the 
Okumura  path  loss  attenuation  relative  to  free  space,  interpolated  from  Figure  10,  pages 
836-840  of  Okumura,  et  al,  for  1775  MHz  with  H  =  45m  and  H  =  3m.  Figure  A-2  gives 

an  approximation  of  figure  A-l,  where  Okumura's  median  attenuation  (in  dB)  relative  to 
free-space  loss  as  a  function  of  distance  is  given  by: 

AFS  =  lloS10d)  (4)/(log102)]  +  33,  ford  <22.  6  km, 
or 

afs  =  [(log1Qd)  (12)/(log102)]  -  3,  for  d  >  22.6  km, 
where 

Apg  =  Attenuation  relative  to  free  space 

4  and  12  (in  dB/octave)  represent  the  slope  of  lines  with  zero  points 
(log  d  =  0;  ie,  d  =  1)  at  33  and  -3  dB,  respectively. 

The  pivot  point,  or  intersection  of  these  two  lines  is  22. 6  km.  This  approximation  is 
almost  identical  for  d  <  15  km  and  d  >  30  km.  For  15  <  d  <30  km,  the  approximation 
is  within  2  dB  of  Okumura's  data. 

However,  the  attenuation  as  given  by  Okumura  is  only  a  median  value.  Thus,  per  Figure 
37  A,  page  860  of  Okumura,  et  al,  a  'correction"  factor  must  l)e  added  to  A  ’  This  is 

l4  u 

accomplished  by  generating  a  random  variate  using  a  normal  distribution  with  a  mean  of  0 
and  standard  deviation  of  5. 

The  total  attenuation,  as  a  function  of  distance  and  frequency  is  as  follows: 

Path  loss  attenuation  (dB)  =  free  space  loss  (dB)  + 

median  Okumura  (dB)  + 
correction  (dB) 


A-l 


DISTANCE 
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A-2.  Approximate  Okumura  Path  Loss  Attenuation  (Relative  to  Free  Space). 


okumura  path  loss  approximation 


Path  loss  attenuation  =  32. 5  +  20  log  F  +  20  log  D 

+  Apg  +  X~n(x,  0,  5) 

The  slopes  of  the  lines  for  figure  A-2  are  independent  of  changes  in  frequency.  However, 
the  slopes  are  dependent  on  the  effective  transmission  height  of  the  repeater.  Even  so,  if 
an  average  effective  transmiss’on  height  can  be  assumed  for  the  repeaters  in  the  packet  radio 
network,  Okumura's  data  can  be  approximated  for  the  N.  A.  C.  simulation  program. 

A  listing  of  the  computer  algorithm  and  sample  input  and  output  data  follows.  The  program 
calculates  the  path  loss  attenuation  for  a  given  distance  and  frequency  based  on  the  above 
loss  equation  and  provides  a  random  component  from  a  normal  distributed  random  number 
generator.  It  is  based  on  antenna  heights  of  repeater  and  terminal  of  45  meters  and  3 
meters,  respectively.  The  only  input  data  that  the  program  requires  is  the  frequency  in 
megahertz  and  distance  in  miles. 
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THIS  PROGRAM  WILL  CALCULATE  THE  PATH  LOSS  ATTENUATION  (RELATIVE  TO  ISOTROPIC)  FOR  A 
TRANSMITTED  SIGNAL  USING  THE  ATTENUATION  DATA  OF  OKUMURA,  ET.  AL  .  ADJUSTED  PER 
OKUMURA  BV  A  NORMALLY  DISTRIBUTED  RANDOM  VARIABLE. 

THE  REPEATER  EFFECTIVE  TRANSMISSION  HEIGHT  IS  ASSUMED  TO  BE  45  METERS  =  14B  FEET  THE 
TERMINAL  EFFECTIVE  TRANSMISSION  HEIGHT  IS  ASSUMED  TO  BE  3  METERS  *  9.B  FEET 


REAL  X(102)  /1 02.0/ 

REAL  DIST,  DISTKM,  FRESP,  OKATEN,  FREQ,  XMEAN,  3TDEV,  TOTATN 
INTEGER  NRAND,  NUMB 

SET  ALL  USER  SUPPLIED  INPUTS  - 

DIST  -  REPEATER  TO  TERMINAL  GROUND  RANGE  (MILES), 
DISTKM  =  SAME  AS  DIST  (KM),  CALCULATED  BY  PROGRAM. 

FREO  *  TRANSMISSION  FREOUE»,CY  (MHZ), 

XMEAN  =  MEAN  OF  RANDOM  NUMBERS  (DB). 

STDEV  =■  STANDARD  DEVIATION  OF  RANDOM  NUMBERS  (DB), 
NRAND  =  NUMBER  OF  RANDOM  NUMBERS  TO  BE  GENERATED, 

X  (1)  =  SEED  NUMBER  FOR  THE  RANDOM  NUMBER  GENERATOR. 

DIST  ■  10.0 
FREQ  -  1775.0 
XMEAN  -  0.0 
STDEV  *  5.0 
NRAND  *  102 
X  (1)  -  (2.0B35)  -  10 

DO  FOR  DIST  -  10  AND  40 

DO  B  J  *  1.2 


CALCULATE  THE  FREE  SPACE  ATTENUATION  (DB) 

FRESP  «  36.6  +  20.0" ALOG10( FREQ  +  20  O'ALOGIOIDIST) 

CALCULATE  THE  OKUMURA  ATTENUATION  (DB)  RELATIVE  TO 
FREE  SPACE  AS  GIVEN  BY  FIGURE  10,  PAGES  B36  40  OF 
OKUMURA,  ET.  AL 

DISTKM  -  (DIST) *(  1 .609344) 

IF  (DISTKM  GT  22  6)  GO  TO  1 

OKATEN  ■  (( ALOG1 0( DISTKM )' (4,0) l/ALOG  10(2,0”  33.0 

GO  TO  2 

1  OKATEN  -  ((ALOG10(DISTKM)'(  12,0) l/ALOG  10(2.0)  3,0 

2  CONTINUE 

IF  (J  GT  1)  GO  TO  3 

GENERATE  A  NORMALLY  DISTRIBUTED  RANDOM  VARIABLE  TO 
CORRECT  THE  OKUMURA  ATTENUATION  PER  FIGURE  37  (A), 
PAGE  B60  OF  OKUMURA,  ET.  AL. 

THE  FIRST  100  NUMBERS  GENERATED  WILL  BE  IGNORED 


K  *  100 
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CALCULATE  THE  TOTAL  ATTENUATION  |DB). 

K  -  K  +  1 

TOTATN  =  FRESP  +  OKATEN  +  X  (K) 

WRITE  THE  TOTAL  ATTENUATION. 

NUMB  *  K  —  100 
IF  (J  GT.  1)  GO  TO  6 
WRITE  (7,4) 

4  FORMAT  (1  HI, //////) 

WRITE  (7,5) 

5  FORMAT  (10X/NUM8ER’,  5X, 'DISTANCE',  5X,'FREE  SPACE',  5X, 

X  'OKUMURA',  5X, 'RANDOM',  5X, 'TOTAL  A ITENUATION',//) 

6  WRITE  (7,7)  NUMB,  DIST,  FRESP,  OKATEN,  X(K),  TOTATN 

7  FORMAT  (12X,  111,  8X,  1F7  2,  7X,  1F7  2,  6X,  1F7.2,  5X,  1F7.2,  10X,  1F7.2,/) 

RESET  DIST 

DIST  =  40.0 

8  CONTINUE 
WRITE  (7,4) 

ALL  DONE. 

STOP 

END 


NUMBER 

DISTANCE 

FREE  SPACE 

OKUMURA 

RANDOM 

TOTAL  ATTENUATION 

1 

10  00 

121  58 

49,03 

.20 

170  42 

2 

40.00 

133  63 

69.10 

2  80 

199.92 

I 

l 
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